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1.0  INTRODUCTION 


Model  displacement  and  the  relative  position  of  the  wings  with 
respect  to  the  fuselage  are  important  measurements  in  aerodynamic 
testing  programs.  The  amount  of  actual  model  pitch,  roll,  and  yaw, 
in  addition  to  magnitude  and  frequency  of  wing  and  fuselage  oscillations, 
are  all  pertinent  to  the  aerodynamic  measurements  of  a test  model 
(Ref.  1).  These  are  especially  critical  in  high  Reynolds  number  test 
conditions  where  aerodynamic  loading  can  be  sufficient  to  produce  large 
static  wing  and  sting  deflections.  Model  deflection  data  must  be  avail- 
able to  test  personnel  both  during  and  following  each  run  to  aid  in  assess- 
ing whether  the  model  or  sting  has  approached  structural  limits  and  to 
provide  the  user  with  information  about  model  deformation  and  position 
in  the  flow  for  corrections  to  the  aerodynamic  data.  During  the  run, 
the  model  may  undergo  programmed  changes  in  the  angle  of  attack  in 
either  a continuous  or  a pitch-pause  mode.  Therefore,  the  position  of 
the  model  relative  to  fixed  reference  or  relative  to  another  model  must 
be  known.  Since  the  measurement  system  must  not  disturb  the  flow, 
conventional  measurement  techniques  are  inappropriate  for  use  in  wind 
tunnels. 

A laser  interferometric  technique  is  proposed  for  these  measure- 
ments based  on  similar  techniques  established  by  others  (Refs.  2 
through  7).  This  technique  requires  only  small  retroreflectors 
mounted  flush  with  the  model  surface  as  preparation  for  meas- 
urement. Resolution  of  the  proposed  technique  in  measuring 
absolute  or  differential  displacements  is  one -quarter  of  the 
wavelength  of  the  laser  illumination  source. 

This  report  describes  the  two-beam  optical  system  and  electronic 
signal  processors  developed  to  implement  the  technique  on  an  experi- 
mental basis.  A description  of  the  two-beam  and  a proposed  multiple- 
beam  optical  system  and  an  analysis  of  the  salient  optical  system 
parameters  are  presented  in  Section  2.  0.  The  electronic  signal  data 
processors  which  provide  the  displacement  information  are  discussed 
in  Section  3.0.  The  errors  inherent  in  this  system  attributable  to 
changing  laser  parameters  and  the  effects  of  variations  in  the  index  of 
refraction  of  the  airflow  are  considered  in  Section  4.0.  Section  5.0 
summarizes  the  work  completed,  and  recommendations  for  future  work 
are  presented. 


Preceding  page  blank 
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2.0  OPTICAL  SYSTEM  ANALYSIS 


An  experimental  two-beam  optical  system  has  been  developed  as  the 
basic  building  block  for  a proposed  multibeam  system.  This  two-beam 
system  provides  operating  criteria  concerning  the  mixing  of  light  from 
retroreflectors,  retroreflector  characteristics,  photodetector  sensi- 
tivity, laser  power  requirements,  and  electronic  signal  processing. 

The  design  of  the  two-beam  system  is  considered  first  followed  by 
the  design  considerations  for  a multibeam  system.  Theoretical  discus- 
sions concerning  design  parameters  which  apply  to  both  systems  are 
then  given. 


2.1  TWO-BEAM  OPTICAL  SYSTEM 

The  optical  system  of  the  two-beam  laser  interferometer  developed 
to  measure  relative  or  absolute  displacements  is  shown  schematically 
in  Fig.  1.  Light  from  a five -milliwatt  helium-neon  laser  is  polarization- 
rotated  90  deg  by  a half-wave  (\/2)  plate  and  reflected  to  a 34 -MHz 
Bragg  cell  by  mirror  M^.  This  rotation  is  necessary  to  allow  the  light 
to  pass  undeflected  through  the  polarization-sensitive  prisms.  In 
passing  through  this  Bragg  cell,  part  of  the  beam  is  up-shifted  in  frequency 
and  the  remaining  part  called  the  zeroth  order  is  unperturbed.  The 
positive  (up-shifted)  first-order  beam  and  the  zeroth  order  are  reflected 
to  a be  am -separating  lens  system  by  mirror  M2.  The  first  order 
beam  is  reflected  by  mirror  M3  and  passes  through  a birefringent 
(Wollaston)  prism  where  it  is  reflected  by  mirrors  M4  and  M5  to  a 
retroreflector  located  on  the  test  surface. 

These  latter  mirrors  (M^  and  M5)  allow  the  separation  and  angle 
between  the  ilium 'nation  beams  to  be  easily  adjusted.  The  unmodulated 
zeroth-order  beam  is  modulated  by  a second  Bragg  cell  which  is 
driven  at  49  MHz  and  oriented  to  enhance  the  positive  first-order 
deflection.  This  Bragg  cell  serves  to  frequency  modulate  a portion 
of  the  light  beam  producing  an  overall  frequency  shift  between  the 
two  beams  of  15  MHz.  The  zeroth  order  or  unperturbed  beam  from 
this  Bragg  cell  is  blocked  to  prevent  extraneous  frequencies  from 
being  generated  should  it  interfere  with  the  other  light  beams. 

Using  two  Bragg  cells,  rather  than  one,  to  obtain  the  desired  15- 
MHz  optical  carrier  facilitates  the  isolation  from  the  optical  carrier 
frequency  of  the  individual  radio  frequency  (RF)  signals  required  to 
drive  the  Bragg  cells  (Ref.  7).  This  permits  locating  the  photodetector 
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and  Bragg  cells  in  close  physical  proximity  without  introducing  cross 
talk  into  the  photodetector  at  the  optical  carrier  frequency. 


Mirror  M7  deflects  the  49-MHz  frequency-shifted  beam  from  the 
second  Bragg  cell  to  a Gian-Air  prism  and  to  a second  retroreflector. 
Total  internal  reflections  within  the  retroreflectors  cause  the  incident 
polarization  vector  to  be  rotated.  The  two  birefringent  prisms 
therefore  separate  the  vertical  polarized  portions  of  the  returned  beams 
and  deflect  them  to  the  beam  splitter.  The  two  beams  are  combined  by 
the  variable  beam  splitter  and  are  reflected  to  the  photodiode  by  mirror 
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Figure  1.  Two-beam  optical  system  schematic 


Mixing  or  optical  heterodyning  the  two  beams  at  the  photodetector 
can  be  thought  of  as  establishing  virtual  fringes,  of  period  \/2,  normal 
to  the  propagation  vector  of  the  retroreflector  illumination  light. 
Because  the  two  beams  are  frequency-shifted  by  15  MHz,  these  fringes 
are  moving  relative  to  a stationary  reference.  Thus  the  light  returned 
by  a stationary  pair  of  retroreflectors  will  be  modulated  at  15  MHz. 
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As  a retroreflector  moves,  the  velocity  of  the  reflected  fringes  will 
appear  to  speed  up  or  slow  down  relative  to  the  direction  of  motion. 
This,  in  turn,  is  manifested  by  a change  in  the  signal  frequency. 

For  clarity,  the  two  retroreflector  illumination  beams  will  be 
referred  to  separately  as  the  object  beam  and  the  reference  beam, 
respectively.  The  output  signal  from  the  photodiode,  a square-law 
detector,  is  proportional  to  the  square  of  the  sum  of  the  incident 
electric  fields,  i.  e.  , 

i(t)  “ {eo  cos  [u0  t + o0<t>]  ■»  ER  cos  [uRt  + oR  (t)]}  2 (1) 

where  uQ  and  uR  are  the  frequencies  >>’  the  object  and  reference 
illuminating  beams  with  amplitudes  i i ( , and  ER,  respectively,  and  <j) 
and  qR  are  the  corresponding  phase  changes  caused  by  optical  path 
length  variations.  Both  reflector  motion,  x(t),  and  the  variations  in 
the  index  of  refraction,  n,  alone  the  optical  path  due  to  density  changes 
in  the  medium  are  included  in  these  phase  terms.  In  general, 

o 7 j n dx  (2) 

where  the  integration  is  along  the  light  path  and  > is  the  wavelength  of 
light  in  a vacuum.  For  a constant  index  of  refraction,  Eq.  (2)  becomes 

0 (t)  = ^ rj  x(t)  (3) 
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Expanding  Eq.  (1)  and  neglecting  the  optical  frequency  terms  to  which 
the  photodiode  is  unable  to  respond,  the  output  signal  in  terms  of  the 
measured  intensity  I is  given  by 

i(t)  a • IQ  + IR  + 2 VlQ  IR  COS  [uct  + Y (%  Xo  (t)  - ^R  XR(t>)](<4> 

where  uc  = ut)  - u is  the  15-MHz  frequency  difference  generated  by 
the  two  Bragg  cells,  I indicates  intensity,  and  the  subscripts,  o and 
R,  indicate  parameters  in  the  object  and  reference  beams,  respective- 
ly. The  first  two  terms  in  Eq.  (4)  establish  a d-c  bias  at  the  output 


the  third  term  rides.  Phase  modulation  is  manifested  on  the 
photodiode  signal  as  a frequency  deviation  given  by 


A u , 
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dq 

dt 


47r  r\  A v(t) 
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(5) 


where  Av(t)  = vu(t)  - vR(t),  vQ(t)  and  vj^(t)  being,  respectively,  the 
components  of  the  object  and  reference  beam  retroreflector  veloci- 
ties in  the  directions  along  the  beam  axes  and  where  r)  is  assumed 
to  be  constant  in  time.  The  effects  of  variation  in  the  index  of 
refraction  generated  by  flow  density  variations  in  the  wind  tunnels 
and  about  the  model  will  be  considered  in  Section  4.2. 


2.2  PROPOSED  MULTIPLE-BEAM  SYSTEM 

An  optical  system  proposed  to  generate  and  mix  multiple  light 
beams  for  more  extensive  and  detailed  interferometric  measurements 
of  model  deformation  and  displacement  is  schematically  shown  in 
Fig.  2.  The  laser  light  is  split  into  two  beams  by  a Bragg  cell  posi- 
tioned to  enhance  the  first  order  and  is  driven  at  49  MHz.  A pair  of 
prisms,  and  P2,  increases  the  separation  between  the  two  beams. 
The  upper  beam  is  down-shifted  by  34  MFz  by  a second  Bragg  cell 
resulting  in  an  overall  frequency  shift  between  the  two  beams  of  15 
MHz.  An  array  of  beam  splitters  and  mirrors  alternately  directs  each 
beam  toward  a series  of  polarization-sensitive  prisms.  (These  could 
be  Wollaston,  as  shown  in  Fig.  2,  or  Gian-Air  prisms.  The  ultimate 
choice  will  depend  on  ease  of  alignment  and  efficiency  in  light  propa- 
gation of  the  optical  system. ) All  but  one  of  the  beams  are  directed  to 
the  model.  One  light  beam  is  directed  toward  a reference  reflector, 
Mp  for  an  absolute  displacement  measurement.  The  rest  are  directed 
to  the  model  through  small  lenses  Lp  of  focal  length,  fp  which  expand 
the  beams  and  direct  them  to  a single  lens,  L2,  with  focal  length  f2. 
Each  small  lens  is  separated  from  L2  by  f^  + f^  which  results  in  a 
collimated  output  beam  with  diameter 


D 

o 


(6) 


where  dp  is  the  beam  diameter  at  Lp  The  small  lenses  are  positioned 
on  a plane  parallel  to  the  model  to  illuminate  each  area  containing  a 
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retroreflector.  Their  separation,  s,  in  this  plane  is  dictated  by  the 
retroreflector  separation,  S,  the  focal  lengths  of  the  lens  system,  and 
the  height,  H,  of  L2  above  the  model,  and  is  given  by 


s g (fj  + f2)s 

H 


(7) 


Light  returned  by  each  reflector  is  collected  and  collimated  by  this 
same  lens  system.  The  resulting  beam  diameter,  dp,  is 


(8) 


where  Dp  is  the  diameter  of  the  effective  aperture  of  the  retroreflector. 
Each  beam  again  enters  the  polarization-sensitive  prisms  and  because 
of  the  polarization  rotation  within  the  retroreflectors  the  perpendicular 
component  is  deflected  toward  the  photodetectors. 


■ Retroreflectors 

Figure  2.  Conceptual  multiple-beam  optical  system. 
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An  array  of  beam -combining  optical  blocks  mixes  adjacent  beams 
and  directs  them  to  the  photodetectors.  To  obtain  equal  intensity 
beams,  the  surfaces  of  these  blocks  are  coated  for  50-percent  trans- 
mission and  reflection.  Each  block  functions  as  two  beam -combining 
elements.  This  is  shown  in  more  detail  in  Fig.  3.  Incident  light  is 
reflected  by  the  front  surface  of  one  block  to  the  back  surface  of  an 
adjacent  block  where  it  is  combined  with  the  light  which  has  passed 
through  that  block.  Block  thickness  is  chosen  to  prevent  secondary 
reflections  from  illuminating  the  photo-sensitive  surface  of  the  detector 
especially  when  the  model  moves  laterally.  Approximately  50  percent 
of  the  light  is  lost  at  the  beam- combining  blocks  due  to  internal  reflec- 
tion and  transmission. 


To  Photodetectors 


Figure  3.  Beam-combining  optical  blocks. 
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Lateral  motion,  L,  of  the  model  will  cause  light  returned  by  the 
reflector  to  move  across  the  face  of  the  photodetector.  The  amount 
of  this  motion,  H,  is  given  by 


(9) 


_3 

As  an  example,  if  f ^ = 0.  2 cm  and  f_  = 50  cm,  then  H - 4 x 10  L. 
Therefore,  if  the  reflectors  were  to  move  laterally  by  5 cm  the  optical 
system  will  reduce  this  motion  at  the  photodetector  to  only  0.  2 mm. 
This  is  an  important  consideration  for  angular  motion  associated  with 
pitch,  roll,  and  angle  of  attack. 

An  array  of  photodetectors  converts  the  optical  information  to 
electrical  signals  which  are  subsequently  processed  for  displacement 
data.  Avalanche  photodiodes  are  proposed  for  this  purpose  because 
of  their  high  sensitivity  (0.2  amps  /watt),  high  gain-bandwidth  product 
(80  GHz),  and  low  system  noise  equivalent  power  (NEP)  of  2 x 10'^ 
watts/  x/llz . In  conjunction  with  a low  noise  amplifier  thece  detectors 
have  a responsivity  of  2 x 10u  volts/ watt  for  a demodulation  bandwidth 
of  50  MHz  (Ref.  8). 


2.3  OPTICAL  SIGNAL-MIXING  REQUIREMENTS 

Optimum  mixing  of  the  light  beams  at  the  beam -combining  blocks 
occurs  when  the  optical  phase  is  uniform  across  each  wavefront.  This 
requirement  is  met  under  the  following  conditions  for  the  two  beams 
(Refs.  7 and  9): 

1.  They  must  have  the  same  mode  structure, 

2.  They  must  be  coincident  and  have  equal  diameters, 

3.  They  must  propagate  in  the  same  direction,  i.  e.  , their 
Poynting  vectors  must  be  coincident  at  the  beam  combiner, 

4.  Their  wavefronts  must  have  the  same  curvature,  and 

5.  They  must  be  identically  polarized. 
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The  first  condition  is  met  by  using  a laser  operating  in  the  funda- 
mental (TEMqo)  mode.  The  second  and  third  condition  can  be  closely 
met  by  fabricating  the  optical  component  mounts  with  the  required 
degrees  of  freedom  for  precise  beam  steering  inside  the  optical  package. 
Focusing  at  the  detector  closely  achieves  similar  wavefront  curvature 
at  the  detector  although  in  a real  optical  configuration,  condition  2 is 
compromised.  The  polarization  prism  ensures  optical  mixing  is 
achieved  with  identical  polarization  vectors.  However,  as  will  be 
seen  in  Section  4.  2,  the  phase  of  the  returned  beams  could  be  dis- 
torted by  density  variations  in  the  flow  about  the  model. 


2.4  RETROREFLECTOR  ANALYSIS 

The  retroreflectors  (corner  cube  prisms),  which  are  mounted 
flush  with  the  surface,  reflect  the  light  back  to  the  source  over  large 
illumination  angles.  Thus,  as  the  surface  moves  due  to  pitching, 
rolling,  or  twisting  of  the  model,  optical  alignment  will  not  be  affected. 

The  intensity,  1^,  of  the  light  that  is  retroretlected,  is  given  by 

W1-'3^)  (10)  I 

where  Ij  is  the  incident  intensity,  i3  is  the  air-glass  interface  loss  due 
to  reflection.  Dp  is  the  diameter  of  the  reflector  aperture,  Dq  is  the 
incident  beam  diameter,  and  A (6^)  is  the  effective  aperture  of  the  corner 
reflector  assuming  perfect  reflecting  surfaces.  Figure  4 illustrates 
the  variation  of  A(0j)  for  a reflector  with  a circular  aperture  and  an 
index  of  refraction  of  1.54  (Ref.  10).  The  maximum  tolerable  angle  of 
incidence,  6^,  will  depend  on  the  minimum  signal  level  required  by  the 
data  processing  electronics.  This  in  turn  is  a function  of  the  sensitivity 
of  the  photodetector  and  output  noise.  Experiments  with  an  avalanche 
photodiode  detector  have  shown  that  can  be  as  much  as  32.  5 deg 
before  the  signal  is  unacceptable  for  the  electronic  signal  processor. 

Total  internal  reflection  within  the  retroreflector  rotates  the 
orthogonal  components  of  the  polarization  vector  which  are  parallel 
and  perpendicular  to  each  reflector's  plane  of  incidence.  Linearly 
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polarized  light  therefore  becomes  elliptically  polarized  on  total  reflec- 
tion. The  amount  of  rotation,  6,  between  the  two  component  vectors 
for  a single  reflection  is  (Ref.  11) 

2 2 
sin“  0.-77 

tan  | = g (H) 

sin“  0. 
l 


Angle  of  Incidence,  0j,  deg 

Figure  4.  Variation  of  the  effective  aperture  of  a circular  retroreflector 
with  angle  of  incidence. 


A maximum  6 of  48  deg  occurs  at  an  incident  angle  of  50.  4 deg  for 
77  = 1.  54.  The  efficiency  of  a polarization-sensitive  optical  system  in 
separating  the  reflected  light  from  the  incident  light  is  therefore  not  as 
large  as  when  a mirrored  surface  is  on  the  object.  Compared  to  a 
first  surface  mirror,  a retroreflector  used  in  conjunction  with  a /.1 4 
plate  and  Wollaston  prism  was  found,  experimentally,  to  be  approxi  - 
mately 40-percent  efficient  in  returning  light  to  the  photodetector. 

The  polarization  rotation  of  the  reflectors  could  be  eliminated  by  coat- 
ing the  reflecting  surfaces  with  silver,  thus  eliminating  the  glass-air 
interface.  However,  now  a >/4  plate  must  be  placed  between  the 
polarization-sensitive  prism  and  mirror  (Fig.  1)  to  rotate  the 
polarization  vector  so  that  the  prism  can  deflect  the  returning  light  to 
the  photodetector. 
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2.5  ANGULAR  MOTION  LIMITATIONS 

The  collimated  beams  illuminating  each  retroreflector  (Fig.  2) 
optimize  the  light  returned  to  the  photodetector  since  areas  of  the 
model  without  a reflector  are  not  illuminated.  However,  the  diameter 
of  each  beam,  D0,  and  the  angle  of  incidence,  0^,  restrict  the  amount 
of  angular  surface  motion,  a,  of  each  reflector.  If  the  beam  diameter 
is  too  small,  or  the  illumination  angle  too  large,  then  a slight  rotation 
can  cause  the  reflector  to  move  out  of  the  illumination  and  the  optical 
signal  will  be  lost.  Four  reflector  rotation  and  illumination  angle 
situations  which  limit  surface  motion  are  depicted  in  Fig.  5.  These 
include  illumination  (1)  normal  to  the  surface,  (2)  in  the  direction  of 
rotation,  (3)  against  the  direction  of  rotation,  and  (4)  in  the  direction 
of  rotation  but  at  a steep  angle.  The  second  illumination  scheme 
maximizes  the  angular  rotation  while  the  fourth  indicates  the  severe 
limits  imposed  when  the  angle  of  incidence  is  too  large.  In  this  case 
the  reflector  moves  completely  out  of  the  opposite  side  of  the  illumi- 
nation beam. 


The  maximum  angular  motion  of  the  reflector  for  cases  2,3,  and 
4 is  given  by 


D 

sin  a = sin  0,  ^cos  0.  - 


1/2 


± cos  0. 
1 


/D 


sm2  0.  + ^ cos  6. 


(12) 


where  0.  is  measured  with  respect  to  the  surface  normal  vector  and  R 
is  the  distance  of  the  retroreflector  from  the  axis  of  rotation,  which 
is  assumed  to  be  perpendicular  to  the  plane  of  the  incidence  angle. 

The  illumination  is  also  assumed  to  be  centered  on  the  reflector  as 
shown  in  Fig.  5.  For  perpendicular  illumination  (0.  = 0),  Eq.  (12) 
reduces  to 


cos  a = 


(13) 


The  second  term  in  Eq.  (12)  is  added  to  or  subtracted  from  the  first 
term  depending  on  whether  motion  is  toward  or  away  from  the  illumi- 
nation as  shown  in  Figs.  5b  and  5c,  respectively.  Maximizing  the 
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rotation  angle  by  increasing  the  illumination  angle  leads  to  the  limita- 
tion shown  in  Fig.  5d.  The  maximum  illumination  angle,  0m,  where 
the  reflector  passes  tangent  to  the  edge  of  the  beam,  is 


COS  0m  = 1 “ ^R 


which  is  seen  to  be  equal  to  the  maximum  rotation  angle,  »m,  for 
normal  illumination.  The  angle  of  rotation  at  the  point  of  tangency  is 


As  the  illumination  angle  becomes  larger  than  this  value,  a is  limited 
by  the  amount  of  retroreflector  travel  in  the  beam  before  it  exits  out 


Retroreflector 

(RR) 


a.  Normal  illumination  (0,  = 0) 


Ij/sL.L. 

k R -RR 

c.  Rotation  opposite  to  the  angle  of 
incidence 


/:>  MM  \ 

H R n''  RR 

b.  Rotation  with  the  angle 
of  incidence 


fa  yi 

— R 


d.  Rotation  through  the  illuminating 
beam 


Figure  5.  Retroreflector  motion  within  a collimated  illumination  beam. 
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of  the  opposite  edge.  The  angle  where  it  initially  leaves  the  beam  is 
given  by 


D 


sin  a = sin  0.-  I cos  9.  + —— 

1 \ i 2R; 


- cos  6. 

l 


9 D /D^ 

sin  9.  - cos  9.  -1^ 
l R l \2R/ 


2 _ 


1/2 


(16) 


Equations  (12)  and  (16)  are  plotted  in  Fig.  6 for  various  values  of 
Dq/R  ranging  from  0.  01  to  0.  5.  Severe  reduction  in  the  rotation 
angle  as  the  illumination  angle  increases  is  quite  evident  in  the  figure. 
As  an  example,  if  a wing  retroreflector  20  in.  from  the  model  center- 
line  is  illuminated  by  a 4-in.  -diam  beam  (DQ/R  = 0.  2)  at  an  incidence 
angle  of  20  deg,  then  the  wing  can  rotate  53  deg  before  the  reflector 
leaves  the  light.  However,  by  increasing  the  incidence  angle  to  30  deg 


Figure  6.  Retroreflector  motion  within  a finite  diameter,  collimated, 
illumination  beam. 

the  reflector  is  restricted  to  a 15-deg  rotation.  It  should  be  noted 
that  angle  rotations  greater  than  32  deg  will  cause  a loss  in  the  optical 
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signal  due  to  the  decreasing  retroreflector  aperture  as  discussed  in 
Section  2.4.  Maximum  pitch  and  roll,  therefore,  will  depend  not  only 
on  the  illumination  beam  parameters  but  also  on  retroreflector  place- 
ment on  the  model  with  respect  to  the  lateral  position  of  the  illumi- 
nation port.  By  centering  the  port  between  the  extreme  retroreflectors' 
positions  on  the  model  (Fig.  2)  the  incident  angle  will  be  minimized. 

Smaller  illumination  beams  and/or  larger  model  rotation  can  be 
used  in  illuminating  the  model  by  positioning  the  light  off-center  on  the 
retroreflector.  However,  motion  of  the  reflector  must  be  known  a 
priori.  This  would  be  especially  advantageous  when  the  model  under- 
goes an  initial  deflection.  DQ  in  Eqs.  (12)  to  (16)  would  then  be  replaced 
by  twice  the  distance  from  the  rellector  to  the  far  edge  of  the  beam. 


2.6  LASER  POWER  AND  PHOTODETECTOR  SENSITIVITY 


The  amount  of  laser  power  necessary  to  illuminate  the  model  and 
be  returned  to  the  photodetector  depends  on  the  minimum  signal-to- 
noise  ratio,  SNR,  required  by  the  signal  processing  electronics.  The 
optical  signal  incident  on  the  detector  can  be  determined  by  considering 
the  losses  incurred  by  the  light  as  it  traverses  each  optical  element. 
This  includes  (1)  the  number,  N,  of  light  beams  generated;  (2)  the 
diameter  of  the  collimated  beams,  Dq;  (3)  the  glass -air  interface 
reflection  loss,  /3,  for  each  uncoated  optical  element  (beam  splitters 
are  coated  to  divide  the  intensity  equally  among  the  N beams  and  the 
beam  splitters  are  coated  for  50-percent  transmission  and  reflec- 
tion, and  all  other  elements  will  be  assumed  to  be  uncoated);  (4)  the 
50-percent  light  loss  at  the  beam-combining  blocks  due  to  secondary 
reflections  and  transmission;  (5)  the  retroreflector  light  reflecting 
efficiency  which  is  given  by  Eq.  (10);  (6)  the  efficiency,  7,  of  the 
polarization-sensitive  prisms  in  deflecting  the  light  to  the  photo- 
detector; (7)  the  efficiency,  e , of  the  second  Bragg  cell,  BC2,  in 
diffracting  the  light  into  the  first  order;  (8)  the  gain,  G,  of  the 
photodetector  at  the  laser  wavelength;  and  (9)  the  minimum  signal  - 
to-noise  ratio  required  by  the  signal  processing  electronics. 

(Noise  on  the  output  signal  is  due  primarily  to  the  inherent  noise, 

Np,  of  the  photodetector  and  the  optical  noise  of  the  light.  The 
latter  noise  source  which  is  a laser  parameter  will  be  neglected.  ) 

The  amplitude  of  the  electrical  signal  out  of  the  photodetector  is 
from  Eq.  (4) 


I 


1 
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where  Ij  and  are  the  intensities  of  each  beam  at  the  photodetector 
which  were  generated  by  Bragg  cell,  BC1.  The  two  beams  are 
assumed  to  have  equal  intensities.  Substituting  all  of  the  above  factors 
into  Eq.  (17)  yields  for  the  output  signal  voltage. 


% “ o- 5 0 * A (8i> -^2 

o 


hi  - /3>  e IT 


where  IT  is  the  initial  laser  intensity. 

L-j 

The  minimum  signal-to-noise  ratio,  SNR,  required  out  of  the 
photodetector  is  found  by  dividing  this  equation  by  the  noise  voltage 
generated  by  the  photodetector  which  is  given  by 

V = (NEP)  G VAf~  (19) 

where  NEP  is  the  noise  equivalent  power  at  the  input  of  the  photodetector 
measured  in  watts-Hz~';  and  Af  is  the  bandwidth  of  the  circuit.  Dividing 
Eq.  (13)  by  (19)  and  solving  for  I yields 

Li 

(SNR)-  N • D 2 -(NEP)  VaF 

j _ Q 

L 0.5  7 (l-fj)9A(0i)DR2V(l-/i)  e 

Table  1 lists  the  laser  power  required  and  the  values  of  the  param- 
eters in  Eq.  (20)  which  were  assumed  for  two  illumination  beam  diameters. 
The  values  of  1^  are  surprisingly  low.  This  is  a result  of  the  high  sensi- 
tivity and  low  noise  of  the  avalanche  photodiode  assumed  for  the  photo- 
detector. Under  actual  operating  conditions,  optical  and  electrical 
noise  can  be  expected  to  increase  the  required  power.  However,  by 
antireflection  coating  the  optical  elements  and  silvering  the  retrore- 
flectors  the  required  laser  power  can  be  minimized. 


3.0  DISPLACEMENT  INTERFEROMETER  SIGNAL  PROCESSORS 


The  output  signal  from  the  photodetector  is  a 15- MHz  carrier  that 
has  been  frequency-modulated  by  the  relative  velocity,  Av(t),  between 
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the  two  retroflectors.  The  signal  frequency  from  Eq.  (4)  is 

f(t)  = f + 2r]Av(t)  ^1) 

\ 

Two  techniques  exist  for  recovering  displacement  information  from 
this  frequency  modulation.  One  is  digital  and  can  best  be  described 
from  a fringe  counting  point  of  view.  This  approach  is  discussed  in 
Section  3.1.  The  second  technique  is  analog  and  involves  converting 
the  signal  frequency  to  proportional  analog  voltage  which,  after  sub- 
tracting a constant  corresponding  to  the  15-MHz  carrier,  can  be  integ- 
rated to  obtain  displacement.  This  type  signal  processor  is  discussed 
in  Section  3.  2.  Section  3.  3 describes  experiments  run  to  compare 
results  obtained  with  both  type  signal  processors. 

Table  1.  Multiple-Beam  Optical  System  Parameters  and 
Required  Laser  Power 

SNR  = 100 


(1-3)  = o.  90 

7 =0.65  (polarization-sensitive  prism  deflection  efficiency) 

A(6  ) = 0.  50  ( 9 . = 20  deg  from  Fig.  4) 

D = 0.  25  in.  (retroreflector  diameter) 

rv 

e = 0.  75  (Bragg  cell  No.  2 diffraction  efficiency) 

NEP  = 5 x 10  w/(Hz)*^  (from  Ref.  8 for  a T1XL74  avalanche 

photodiode  module) 
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3.1  PHASE  QUADRATURE  DIGITAL  DATA  PROCESSOR 


This  approach  to  processing  the  photodetector  signal  uses  a digital 
up/down  counter  to  measure  surface  position  by  counting  digital  pulses 
developed  for  each  plus  or  minus  A/4  of  surface  motion.  Knowledge  of 
when  the  surface  changes  direction  is  needed  to  accurately  track  sur- 
face position  by  directing  A/4  pulses  to  the  appropriate  up  oi  down 
counter  terminal.  A pulse  separation  circuit  which  utilizes  the  phase 
difference  between  the  Bragg  cell  and  photodiode  signals  to  instanta- 
neously separate  up  and  down  position  pulses  will  be  described. 

Figure  7 is  a block  diagram  of  the  circuit  used  to  digitally  process 
the  photodetector  signal.  Not  shown  are  the  34-  and  49-MHz  crystal 
oscillators  used  to  drive  the  Bragg  cells.  The  15-MHz  frequency- 
modulated  (FM)  signal,  S(t),  from  the  photodetector,  and  the  34-MHz, 
dj(t),  and  40  MHz,  d2(t),  signals  which  drive  the  Bragg  cells  to  pro- 
duce the  15-MHz  frequency  difference  between  t lie  object  and  the 
reference  light  beams  are  the  circuit  inputs.  These  can  be  repre- 
sented by 


SCt)  = VQ(t)  cos  [«l5t  + ~ Ax(t)J 
S [U34l] 


dj(t) 


V i cos 


(22) 


,(t)  - V2  cos  [«4Bt] 


where  the  subscripts  on  u indicate  the  frequency  value  in  megahertz  and 
assuming  n = 1.  The  amplitude  of  S(t)  is  shown  as  a function  of  time  due 
to  the  change  in  incident  angle  and  subsequently  a change  in  the  returned 
light  amplitude  from  the  retroreflector  as  the  model  moves.  Limiter, 
A2,  provides  a constant  signal  amplitude  for  V The  34-  and  49-MHz 
signals  are  multiplied  together  by  a double -balanced  mixer  to  give  sum 
and  difference  frequencies,  i.  e.  , 


V1V2  V1V2 
g(t)  - — ^ COS  (li-j-t)  + — COS  (i. 


83 


(23) 


A 20 -MHz  lowpass  filter  removes  the  high  frequency  component  and 
an  amplifier,  A3,  separates  the  15-MHz  reference  signal  into  two  com- 
plimentary signals  180  deg  out  of  phase.  The  phase  of  the  noninverted 
signal  is  decreased  90  deg  by  a delay  line  so  that  the  total  phase  shift 
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between  the  two  is  90  deg.  These  signals  can  be  represented  as 

V V 
12 

rl<t)  = — g cos  (^15t  - 7T) 

(24) 

V V 

r2(t)  = g 2 COS  ^15*  ' 

Two  double -balanced  mixers  multiply  these  signals  with  the  photodetector 
signal,  S(t),  to  produce  Mi<t)  and  M2(t),  respectively,  which  after 
algebraic  and  trigonometric  simplification,  can  be  expressed  as 

Ml(t)  ^ -^p~^[cos  ^>3ot  +-y  Ax(t))  + cos  Ax(t)^J 


(25) 


(*)  " "41  2 sin  (u3Qt  + ' sin  (y  Ax(t)) ] 


Lowpass  filtering  the  output  from  each  of  the  mixers  removes  the  30 -MHz 
frequency  term  and  yields 


- V V V / \ 

Pl(t)  = -°  £ COS  Ax(t)J 

-V0V1V2  Uw  \ 

P2(t)  = 4 sin  \X  Ax(t)/ 

Note  that  the  delays  in  the  filters  are  matched  since  subsequent  signal 
processing  will  depend  on  the  phase  between  these  two  signals.  Also, 
note  that  the  5-MHz  break  frequency  of  these  filters  limits  the  maximum 
allowable  surface  velocity  to  ±1.  6 meters/sec. 

It  can  be  seen  from  Eq.  (26)  that  P^(t)  has  the  same  sign  regardless 
of  the  sign  of  Ax(t),  which  changes  with  the  surface  direction.  However, 
when  the  sign  of  Ax(t)  changes,  the  sign  of  P2(t)  also  changes  and  it  is 
this  signal  in  conjunction  with  P^(t)  that  is  used  to  direct  position  pulses 
to  the  up  or  down  terminal  of  the  counter. 

Signals  P^  and  P2  are  converted  to  Transistor  Transistor  Logic 
(TTL)  pulse  trains,  A and  C,  respectively,  in  Fig.  7 by  zero-crossing 
detectors  A4  and  A5.  Delaying  each  signal  by  15  and  30  nsec  yields  TTL 
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(26  a) 


(26  b) 


Figure  7.  Phase  quadrature  position  data  processor. 


pulse  trains  and  B,  respectively.  The  AND  gates,  to  G4,  and  OR 
gates,  G5  and  G6,  implement  the  Boolean  logic  expressions 

U = A • B • C6+  A-B.C6  (27a) 

D = A • B • C + A-B-C  (27b> 

6 6 

where  A,  B,  and  are  the  Boolean  complements  of  the  A,  B,  and 
pulse  trains,  respectively.  These  expressions  generate  position  pulses 
on  either  the  up  or  down  line  (U  or  D expression,  respectively)  of  the 
counters  depending  upon  the  direction  of  the  surface  motion.  Each  posi- 
tion pulse  is  30  nsec  wide  and  is  produced  for  each  X/4  of  surface  move- 
ment. 

The  operation  of  this  logic  circuit  (and  the  Boolean  expressions)  in 
detecting  the  phase  between  the  two  signals  can  be  understood  by  analyz- 
ing the  time  sequence  of  the  waveforms  shown  in  Fig.  8.  The  Boolean 
logic  identity  A • B • C = (A  • B)  • C will  be  utilized  to  show  how  the 
30-nsec  pulsewidths  are  generated.  It  is  seen  in  Fig.  8 that  the  logical 
multiplications  (AND  operations)  of  A with  B and  of  A with  B result  in 
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Figure  8.  Analog  waveforms  and  logic  pulse  trains  for  the  phase  quadrature 
signal  processor. 

30-nsec  pulses  which  occur  at  the  positive  and  negative  going  zero  cross- 
ings, respectively,  of  the  Pi  waveform.  These  position  pulses  are  then 
separated  according  to  the  logic  state  of  C since,  according  to  Eq.  (26b), 
the  sign  of  P2  and  therefore  the  logic  state  of  this  signal  is  dependent  on 
the  sign  of  Ax(t).  (Note  that  the  C pulse  train  is  delayed  by  15  nsec  to 
best  center  and  the  30-nsec  position  pulses  in  order  to  optimize  fre- 
quency response. ) 

To  fully  evaluate  the  pulse  separation  circuit  just  described,  it  is 
necessary  to  verify  that  no  pulses  are  lost  or  gained  on  the  U or  D logic 
line  at  surface  turnaround  for  each  possible  phase  relationship  between 
the  A and  C waveforms.  As  explained  earlier,  the  surface  can  be 
assumed  to  be  cutting  virtual  fringes  with  Pi  and  P2,  the  representative 
expressions  for  the  resulting  signals.  Therefore,  a physical  understand- 
ing of  the  effects  of  surface  motion  and,  correspondingly,  the  phase  dif- 
ference between  Pi  and  P2  on  the  position  pulses  preceding  and  follow- 
ing a direction  change  at  various  locations  in  this  fringe  pattern  can  be 
obtained.  Figure  9 illustrates  zero  crossings  of  Pi  (solid  lines)  and  P2 
(dashed  lines)  along  the  light  propagation  path,  z.  Paths  1 through  4 
indicate  four  possible  phase  relationships  between  the  logic  signals  A 
and  C at  the  turnaround  point.  These  signals  and  the  resulting  pulse 


26 


* 


AEDC  TR-76  116 


Zero  Crossing  of  Pj 
Zero  Crossing  of  P2 
Logic  1 State  of  A 
Logic  1 State  of  C 


Paths  of 
Surface 
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Figure  9.  Pictorial  representation  of  surface  motion-generating  phase-delayed 
signals  P1  and  P2. 

trains  on  the  up  and  down  counting  busses  for  paths  1,  2,  and  3 are 
illustrated  in  Fig.  10.  Signals  for  the  fourth  path  correspond  to  those 
presented  in  Fig.  8 as  a function  of  time.  In  Fig.  10  the  signals  have 
been  shifted  for  coincident  turnaround  points.  A symmetrical  count  is 
observed  in  all  instances,  i.  e.  , there  are  as  many  pulses  on  the  up 
logic  line  as  there  are  on  the  down  line  for  equidistant  motion. 

Noise  on  the  input  signal  S(t)  can  affect  the  count  through  the  A and 
C waveforms  by  adding  or  dropping  pulses  from  the  waveforms.  Pulse 
addition  could  result  in  erroneous  counts  being  accumulated.  However, 
the  probability  of  pulse  addition  or  cancellation  occurring  with  this  cir- 
cuit is  small.  For  example,  if  a noise  pulse  occurred  in  the  A waveform 
while  the  C waveform  remained  either  a logic  1 or  a logic  0,  one  up  and 
one  down  count  pulse  would  be  generated.  The  overall  effect  of  the  noise 
pulse  would  therefore  be  zero.  A A/4  error  would  occur  in  tJ.e  data  only 
if  the  contents  of  the  counters  were  recorded  between  the  rising  and  fall- 
ing edges  of  the  added  pulse.  A similar  argument  holds  for  pulse  cancel- 
lation. Only  when  the  addition  or  cancellation  of  a pulse  by  noise  occurs 
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Figure  10.  Logic  waveforms  of  surface  motion  depicted  in  Fig.  9. 

on  A at  the  same  time  as  a state  transition  on  C or  vice  versa  will  erron- 
eous counts  be  unequally  accumulated.  The  amount  of  this  error  would 
then  be  a function  of  the  input  signal -to-noise  ratio.  Noise  pulses  on  the 
C waveform  could  cause  a position  pulse  to  be  on  the  wrong  counter  input 
only  if  it  occurred  at  the  same  time  as  the  position  pulse  and  was  at 
least  30  nsec  wide. 

Resolution  of  the  displacement  data  can  be  reduced  from  X/4  to  X/2 
thus  reducing  the  number  of  logic  circuit  components  required  and 
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doubling  the  maximum  surface  velocity  capability  of  the  logic  circuitry. 
This  is  accomplished  by  implementing  only  one  term  from  each  of  the 
expressions  in  Eq.  (2  7).  Of  the  four  possible  combinations  for  the  U and 
D functions,  only  two 


u = 

A • 

B • 

c«- 

D 

= A 

■ b c4 

(28) 

u = 

A- 

B • 

D 

= A • 

B'ea 

(29) 

are  appropriate.  The  other  two  pairs  have  the  possibility  of  a pulse 
being  dropped  each  time  the  surface  changes  direction.  As  an  example, 
referring  to  the  waveforms  in  Fig.  8,  and  assuming  that  the  logic  ex- 
pressions U - A • B • Cg  and  D = A . B . are  to  be  implemented,  the 
phase  relationship  of  signals  A and  C is  such  that  a position  pulse  is  not 
present  on  the  up  input  since,  prior  to  surface  turnaround,  the  A wave- 
form is  in  a high  logic  state.  The  expressions  of  Eqs.  (28)  and  (29)  on 
the  other  hand  are  consistent  in  generating  both  an  up  and  a down  pulse 
when  the  surface  changes  direction. 

A pulse  separation  circuit  implementing  both  terms  of  Eq.  (27)  to 
obtain  A/4  resolution  was  fabricated.  Presented  in  Fig.  11  are  photo- 
graphs of  circuit  waveforms  taken  from  an  oscilloscope  to  verify  proper 
circuit  operation.  Figure  11a  shows  Pj(t)  and  P2<t)  at  a surface  turn- 
around point  (horizontal  screen  center)  along  with  their  corresponding 
TTL  logic  waveforms,  A and  C.  Figure  lib  illustrates  formation  of 
the  term  A ■ B • C^from  Eq.  (27).  Figure  11c  illustrates  the  transfer 
of  position  pulses  between  the  up  and  down  counter  lines  at  surface 
turnaround. 

It  is  interesting  to  note  that  this  data  processing  technique  is 
similar  in  operating  principle  to  the  optical  quadrature  technique  which 
utilizes  two  photodetectors  optically  separated  by  90  deg  on  the  inter- 
ference pattern  (Ref.  12).  However,  the  electrical  circuit  requires  half 
as  many  photodetectors  and  is  not  sensitive  to  photodetector  position. 

3.2  ANALOG  DATA  PROCESSOR 

An  analog  approach  used  to  determine  position  from  the  frequency - 
modulated  photodetector  signal  involves  converting  the  signal  frequency 
to  a proportional  analog  voltage,  removing  from  it  a constant  voltage 
proportional  to  the  carrier  frequency,  fc,  and  integrating.  The  problem 
with  this  approach  occurs  in  trying  to  obtain  frequency-to-voltage  con- 
version at  15  MHz  while  still  maintaining  the  accuracies  required.  A 
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solution  to  this  problem  utilizes  frequency  scaling  to  move  the  photo- 
detector signal  to  a lower  frequency  range  which  is  compatible  with 
standard  and  very  accurate  frequency-to-voltage  (f/v)  converters. 

The  converter  used  in  the  prototype  analog  processor  decodes  fre- 
quencies ranging  from  10  Hz  to  110  kHz  to  voltages  between  1 mv  and 
11  v,  respectively.  Typical  linearity  is  ±0.006  percent  with  a maximum 
tolerance  of  ±0.05  percent.  The  converter  time  constant  is  0.05  msec 
which  results  in  less  than  a 0.  008  percent  amplitude  variation  of  output 
voltage  at  frequencies  between  0 and  100  Hz. 


c.  Transfer  of  position  pulses 
at  surface  turnaround 


Figure  11.  Circuit  waveforms  from  the  X/4  resolution  digital  processor. 

Figure  12  is  a block  diagram  of  the  analog  position  measurement 
system.  The  amplitude -limited  (by  limiter  A2),  bandpass -filtered  FM 
signal  is  frequency-translated  to  55  kHz,  the  center  frequency  of  the 
f/v  converter.  It  is  simultaneously  scaled  so  that  the  range  of  frequency 
variations,  Af,  does  not  exceed  the  full-scale  range  of  the  converter. 
Frequency  translation  is  achieved  by  heterodyning  this  signal  with  a 
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Figure  12.  Analog  readout,  frequency-to-voltage  converter. 


higher  frequency  crystal -controlled  oscillator  signal  of  frequency  fn. 

The  sum  frequency  component  is  removed  from  the  difference  frequency 
by  a lowpass  filter.  An  oscillator  higher  in  frequency  than  15  MHz  is 
used  because  it  gives  a greater  separation  between  the  sum  and  differ- 
ence frequencies  than  would  a lower  frequency  oscillator,  thus  easing 
the  rolloff  requirements  on  a lowpass  filter.  However,  the  higher  fre- 
quency oscillator  produces  an  inversion  in  the  direction  of  frequency 
deviation  of  the  translated  FM  signal..  The  output  from  the  lowpass 
filter  is  converted  to  a TTL  compatible  signal  by  zero-crossing  detector 
A3.  Frequency  scaling  is  then  performed  by  digital  division  with  a 
chain  of  flip-flops. 

The  amount  of  frequency  translation  and  division  can  be  optimized 
to  the  full  range  of  the  f/v  converter  with  respect  to  the  maximum  sur- 
face velocity  expected.  This  allows  the  sensitivity  of  the  technique  to 
be  maximized.  However,  55  kHz  minus  the  sum  of  the  highest  modulat- 
ing frequency  (5  kHz)  and  the  maximum  value  of  Af  must  always  be 
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positive  if  the  displacement  information  is  to  be  recovered  without  dis- 
tortion. This  requirement  can  be  met  here  by  limiting  the  maximum 
allowed  value  of  Af  to  50  kHz. 


Frequency  translation  and  division  with  a resultant  center  frequency 
of  55  kHz  can  be  represented  by  the  equation 


(f  - 15  MHz)/2n  = 55  kHz 


n 


(30) 


where  n is  the  number  of  flip-flops  in  the  frequency  divider  and  fn  is  the 
frequency  of  the  translation  oscillator.  The  requirement  that  Af  be 
limited  to  50  kHz  following  frequency  division  dictates  (for  the  maximum 
frequency  deviation)  that 

50  kHz  (31) 


Af  =2 
max 


n 


Since  Afmax  is  a function  of  peak  surface  velocity,  vmax,  given  by 


Af  = (2 /X)  v 
max  max 


(32) 


Equation  (31)  can  be  solved  for  n in  terms  of  vmax  and  X yielding 


n 


where 


[] 


J-  log  2 . 


(33) 


indicates  that  n is  the  next  greater  integer  of  the  value  in 


these  symbols.  After  determining  n,  f is  obtained  from  Eq.  (30)  where 


f = (2 
n 


n 


55  kHz)  + 15  MHz 


(34) 


As  an  example,  assume  that  the  maximum  expected  velocity  is 
1 m/sec.  Substituting  into  Eq.  (33)  with  X = 6328  A gives 


« = [[5-  981]  ‘ 


Using  this  value  of  n in  Eq.  (34)  gives  a value  for  fn  of  18.  520  MHz.  The 
frequency  deviation  from  15  MHz  for  a velocity  of  1 m/sec  is,  from 
Eq.  (32),  3.  16  MHz.  This  means  that  the  largest  desired  frequency  out 
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of  the  balanced  mixer  is  the  difference  between  18.  52  MIIz  and  11.  84  MHz 
(15  MHz  minus  3.  16  MHz)  or  6.  68  MHz.  The  break  frequency  of  the  low- 
pass  filter  should  then  be  greater  than  6.  68  MHz  with  high  attenuation  near 
30  MHz  in  order  to  reject  the  sum  frequency  from  the  balanced  mixer  output. 

Additional  ranges  with  full-scale  velocities  of  0.  5 m/sec,  0.  25  m/sec, 
etc.  can  also  be  made  available  as  shown  in  Fig.  12  by  switching  in  differ- 
ent frequency  translation  oscillators  (16.  76  MHz,  15.  88  MHz,  etc.)  and 
correspondingly  different  values  of  n (5,  4,  etc).  Maximum  sensitivity 
occurs  on  the  range  where  n = 0 and  fn  = 15.  055  MHz.  Here  a velocity  of 
only  ±1.  75  cm/ sec  produces  a full-scale  output  of  10  v from  the  frequency  - 
to-voltage  converter.  Conversely,  values  of  n larger  than  six  with  a 
higher  value  for  the  filter  break  frequency  and  with  higher  frequencies 
for  fn  would  permit  full-scale  velocities  greater  than  1 m/sec  if 
required. 

The  output  of  the  frequency-to-voltage  converter  is  lowpass -filtered 
to  reduce  noise  on  the  velocity -analog  signal.  For  the  type  of  FM  demodu- 
lation system  used,  the  total  system  noise  bandwidth  is  that  of  this 
filter  (assuming  acceptable  input  FM  signal-to-noise  ratios)  (Ref.  13). 

This  filter  also  reduces  ripple  at  the  converter  input  frequency  which 
appears  on  the  output  and,  although  in  the  millivolt  range,  can  contribute 
an  error  into  the  integrated  data. 

Care  must  be  taken  to  ensure  that  the  filter  response  does  not  roll 
otf  significantly,  and,  in  turn,  vary  the  amplitude  of  analog  voltages 
having  frequency  components  below  the  maximum  expected  displacement 
frequency  component  (typically  100  Hz).  The  rolloff  effects  of  various 
practical  filters  were  analyzed  with  a computer.  An  inductor -capacitor 
(LC)  pi  filter  with  a break  frequency  of  approximately  1 kHz  and  with  a 
l-kS7  termination  impedance  was  found  to  be  the  best  of  those  studied. 

A practical  0.3-h  inductor  having  35  ohms  of  series  resistance,  and  a 
pair  of  0.  1 5 -/u f capacitors  in  this  configuration  were  found  to  vary  the 
amplitude  by  only  0.015  percent  below  100  Hz  while  providing  more  than 
40  db  of  attenuation  above  5 kHz.  Insertion  loss  in  the  filter  was  calcu- 
lated to  be  1.  72  percent. 

Operational  amplifier  (op-amp)  A4  (Fig.  12)  removes  the  d-c  term 
corresponding  to  the  55 -kHz  translated  value  of  the  15-MHz  carrier  and 
provides  inversion  so  that  the  polarity  of  the  voltage  into  the  integrator 
corresponds  to  the  direction  of  motion.  Additionally,  since  the  analog 
integrator  will  normally  accept  full-scale  inputs  different  from  the  full- 
scale  output  of  the  converter  and  because  of  signal  attenuation  in  the 
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filter  termination  network,  the  amplifier  gain  is  adjusted  to  optimize  and 
calibrate  the  analog  circuit. 


Calibration  of  the  analog  measurement  system  for  values  of  n and  fn 
is  simply  a matter  of  adjusting  the  gain  of  op-amp  A4  to  a convenient 
value.  This  can  be  done  with  a frequency  counter,  signal  generator,  and 
voltmeter.  If  the  filter  and  op-amp  are  considered  to  have  a composite 
gain  of  G',  the  governing  calibration  equation  is 


Av(t)  = 


V (t) 
v 


(35) 


for  the  velocity  signal  where  Vv(t)  is  the  voltage  out  of  op-amp  A4,  and 


Ax(t)  = 


Gj  G' 


V (t) 
x 


(36) 


for  the  displacement  signal.  Vx(t)  is  the  voltage  at  the  output  of  the  inte- 
grator and  Gj  is  the  integrator  gain. 


A prototype  analog  data  processor  was  fabricated  with  five  full-scale 
velocity  ranges  of  from  ±.l.  6 cm/ sec  to  ±1  m/sec,  each  of  which  pro- 
vides a ±10-v  output  at  full  scale.  To  verify  proper  operation  of  the 
processor,  the  15 -MHz -centered  frequency-to-voltage  conversion  char- 
acteristic of  each  range  was  checked  with  a signal  generator,  frequency 
counter,  and  digital  voltmeter.  In  each  case  the  characteristic  proved 
to  be  a straight  line  within  the  measurement  accuracy  of  the  3.  5-digiL 
voltmeter  used.  The  frequency  response  of  the  output  lowpass  filter  was 
also  checked  and  found  to  agree  with  the  desired  theoretical  response. 


3.3  EVALUATION  OF  DIGITAL  AND  ANALOG  PROCESSORS 

An  experiment  was  conducted  comparing  the  phase  quadrature  digital 
readout  unit  to  a known  standard.  A 4-in.  -diam  micrometer  with  a 1-in. 
travel  was  used  to  push  a gauge  block  having  a retroreflector  attached. 
The  results  were  as  expected  with  the  digital  readout  giving  the  correct 
value  for  the  micrometer  reading  to  a better  resolution  than  could  be 
obtained  from  the  micrometer  scale. 

The  analog  readout  was  calibrated  against  a known  frequency  devia- 
tion and  a frequency-to-voltage  factor  established  as  discussed  in  the 
previous  section.  To  further  evaluate  the  two-beam  system,  a compari- 
son between  the  two  different  readout  techniques  was  undertaken. 
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To  experimentally  compare  the  performance  of  the  digital  and  analog 
data  processors,  displacement  measurements  were  taken  of  a vibrating 
bar  (Fig.  13)  as  a simulation  of  model  motion  in  a wind  tunnel.  The  bar 
was  mounted  horizontally  on  a small  optical  bench  such  that  its  length 
could  be  varied  to  change  the  natural  frequency  of  oscillation.  An  initial 
displacement  was  set  with  a standoff  holding  the  bar  from  the  neutral 
position.  Varying  the  length  of  the  standoff  gave  free  choice  of  the  initial 
displacement.  Rapid  removal  of  the  standoff  set  the  bar  oscillating  about 
the  neutral  position.  Judicious  addition  of  weight  to  the  bar  allowed  some 
control  of  the  oscillating  frequency  and  its  harmonic  modes.  Near  the 
free  end,  a 0.  5-in.  retroreflector  was  mounted  with  the  option  of  mount- 
ing the  other  retroreflector  on  the  oscillating  bar  or  fixing  it  to  a station- 
ary reference. 
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Figure  13.  Experimental  setup  for  the  vibrating  bar  experiment. 


Data  from  the  two  processors  were  recorded  simultaneously  with  a 
minicomputer  data  acquisition  system.  Pairs  of  readings  were  taken  at 
fixed  intervals  of  time  as  determined  by  a gated  pulse  generator.  A flow 
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chart  is  shown  in  Fig.  14  describing  the  data  acquisition.  At  the  beginning 
of  each  data  run,  the  computer  would  await  a command  pulse  from  the  pulse 
generator  to  record  the  first  pair  of  data  points.  The  pulse  generator  was 
gated  ON  just  prior  to  removing  the  standoff  holding  the  deflected  bar.  The 
computer  then  recorded  data  in  pairs  each  time  a pulse  was  generated  until 
the  required  number  of  points  for  the  run  had  been  obtained. 


Figure  14.  Flow  chart  for  data  processing  of  vibrating  bar  experiment. 


I 


Data  from  the  digital  processor  up/down  counter  were  strobed  into 
a storage  register  and  made  available  to  a computer  parallel  data  input 
card  each  time  a pulse  indicated  that  a data  pair  was  to  be  recorded. 
Simultaneously,  an  analog  sample/hold  circuit  tracking  the  analog 
processor  output  was  switched  to  hold  and  conversion  was  initiated  in  a 
14-bit  analog-to-digital  converter.  After  a short  delay  (100  /usee)  to 
allow  the  digital  data  lines  to  settle,  a computer  flag  was  set  on 
the  digital  processor  input  channel.  Similarly,  as  soon  as  conver- 
sion of  the  analog  reading  was  completed,  a flag  was  set  on  the 
analog  data  channel.  When  both  flags  were  set,  the  data  pair  was 
stored  as  equally  indexed  elements  of  two  memory  arrays  reserved 
for  the  respective  processors.  Since  storage  of  the  data  pair  re- 
quired less  than  30  /usee,  the  maximum  rate  at  which  data  pairs 
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could  be  taken  (i.  e. , the  maximum  frequency  to  which  the  pulse 
generator  could  be  set)  was  essentially  determined  by  the  digital 
data  line  setting  delay  of  100  jusec.  Sample  rates  of  up  to  10  kHz 
were  therefore  possible  although  most  of  the  data  taken  were 
recorded  at  rates  of  500  Hz  or  less. 


Memory  arrays  for  the  data  were  reserved  by  a supervisory  BASIC 
language  program.  After  setting  up  the  array  allocation,  the  BASIC 
program  transferred  control  to  an  assembly  level  subroutine.  This 
subroutine  acquired  the  data  in  pairs  each  time  the  two  flags  were  set. 
Once  the  pulse  generator  was  gated  ON  and  the  number  of  pairs  request- 
ed by  the  BASIC  program  was  obtained,  the  assembly  routine  converted 
the  data  to  the  format  required  by  the  BASIC  compiler.  The  arrays  and 
control  were  then  returned  to  the  supervisory  BASIC  program  for  all 
further  operations. 

BASIC  level  programming  was  used  to  store  each  data  array  on  mag- 
netic tape  for  future  reference  and  to  convert  the  digital  array  from 
counts  to  meters  and  the  analog  array  from  volts  to  meters  per  second. 
Similar  programming  was  used  to  integrate  the  analog  velocity  data  and 
to  scale  and  plot  the  data  via  an  online  graphics  terminal  thus  providing 
the  data  curves  given  here. 

For  the  data  presented,  900  simultaneous  analog  (velocity)  and 
digital  (displacement)  points  were  taken  and  plotted.  A trapezoidal 
approximation  of  the  area  under  every  four  points  of  the  velocity  data 
was  used  to  obtain  225  integrated  velocity  (analog  displacement)  points. 

It  should  be  noted  that  since  fewer  integrated  velocity  points  were  avail- 
able for  plotting,  a somewhat  coarser  curve  results.  Also,  some  error 
should  be  expected  in  the  integration  process  when  high  frequencies  are 
contained  in  the  velocity  data. 

A large  variety  of  bar  oscillation  runs  were  taken.  Only  three  runs 
will  be  presented  to  show  the  dynamic  range  of  the  interferometer  in 
addition  to  presenting  the  data  from  two  different  readout  techniques: 
one  from  the  digital  processor  (fringe  counting),  and  the  other  from  the 
analog  processor  (frequency-to-voltage  conversion  of  the  photodetector 
signal).  The  first  two  experiments  that  will  be  discussed  use  the  pri- 
mary optics  module  only,  without  beam  expanding.  The  last  run  to  be 
discussed  used  a beam-expanding  module  to  produce  1-in.  beams. 

The  first  run  is  shewn  in  Fig.  15.  Both  retroreflectors  were 
mounted  7 in.  apart  on  the  oscillating  bar,  and  initial  displacement  of 
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the  bar  was  set  at  approximately  1.  65  cm.  The  bar  was  strongly  damped 
to  reduce  the  fundamental  and  harmonic  modes.  Figure  15a  displays  the 
relative  velocity  between  the  two  retroreflectors.  This  figure  also  shows 
some  residual  "background"  vibrations  of  the  bar  from  acoustic  excita- 
tion. Computer  integration  yields  Fig.  15b  which  is  relative  displace- 
ment between  the  two  retroreflectors.  The  initial  peak  excursion  is 
seen  to  be  approximately  510  p with  an  immediate  exponential  decay  to  a 
final  position  of  approximately  31  p.  Direct  fringe  counting  yields  a dis- 
placement in  steps  of  X/4  or  0.  16  p.  Figure  15c  is  the  plot  of  the  fringe 
counting  output  (digital)  which  may  now  be  compared  to  the  integrated 
velocity  data  of  Fig.  15b.  Excellent  agreement  is  seen  between  the  two 
readout  methods  with  almost  identical  peak  amplitude  excursions  even 
though  a "rough"  integration  technique  is  used  in  the  analog  method. 
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a.  Online  velocity  data  from  analog  readout 
Figure  15.  Data  from  Run  No.  1. 
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A larger  displacement  is  seen  in  the  second  run  (Fig.  16)  with  strong 
damping  of  the  fundamental  and  moderate  damping  of  the  harmonic  fre- 
quency. The  same  experimental  setup  was  used  as  in  the  first  run.  The 
two  frequencies  are  easily  deduced  from  the  velocity  data.  Fig.  16a,  of 
the  analog  processor.  The  online  computer  integration  is  shown  in  Fig. 
16b.  Initial  peak  excursion  is  seen  to  be  0.  282  cm  with  only  a minor 
contribution  from  the  higher  frequency  component.  Figure  16c  shows  the 
displacement  from  the  direct  fringe  counting  processor  and  is  seen  to  be 
virtually  identical  to  the  integrated  velocity  data.  The  fringe  counting 
displacement  also  shows  a peak  amplitude  of  0.  282  cm  identical  to  the 
integrated  velocity.  Since  the  fringe  counting  displacement  data  have  a 
0.  16 -/a  resolution,  the  data  represented  in  the  last  25  percent  of  the 
data  rim  may  be  retrieved  and  analyzed.  Figure  16d  shows  the  last  25 
percent  of  Fig.  16c  easily  showing  the  fundamental  frequency  with  approxi- 
mately 35 -/a  amplitude  and  the  harmonic  at  approximately  5-jU  amplitude 
thereby  demonstrating  the  resolution  available  in  large  amplitude  data. 


Time,  sec 

a.  Online  velocity  data  from  analog  readout 
Figure  16.  Data  from  Run  No.  2. 
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b.  Integrated  velocity  of  Fig.  16a 


c.  Online  displacement  from  digital  readout 
Figure  16.  Continued. 
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d.  Instant  recall  of  last  25  percent  of  Fig.  16c. 

Figure  16.  Concluded. 

To  demonstrate  the  large  amplitude  and  velocity  excursion  capabil- 
ity of  these  readouts,  the  third  run  was  made  with  the  oscillating  bar 
extended  to  its  maximum  length  of  60  in.  and  with  one  retroreflector  ref- 
erenced to  a stationary  surface.  This  setup  allowed  for  direct  displace- 
ment measurement  of  the  tip  of  the  oscillating  bar.  A highly  damped 
fundamental  was  obtained  with  minimal  damping  of  the  harmonic  mode  as 
shown  in  Fig.  17.  A beam -expanding  output  module  was  incorporated  in 
this  experiment  because  of  the  lateral  motion  of  the  oscillating  retro- 
reflector.  Peak  excursion  of  the  bar  was  measured  to  be  20  cm.  Again, 
excellent  agreement  between  the  two  displacement  measurements  is 
seen  demonstrating  the  dynamic  range  of  the  interferometer  and  applic- 
ability to  a wind  tunnel  test. 
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c.  Online  displacement  from  digital  readout 
Figure  17.  Concluded. 


4.0  DISPLACEMENT  ERROR  ANALYSIS 


4.1  LASER  FREQUENCY  AND  OPTICAL  SYSTEM  STABILITY 


The  sensitivity  of  the  interferometer  for  displacement  measurements 
also  makes  it  susceptible  to  motion  of  the  optical  system  components  and 
to  frequency  variations  in  the  output  of  the  laser.  A laboratory 
interferometer  operating  on  a stable  table  in  a controlled  environ- 
ment has  a theoretical  accuracy  dependent  on  the  stability  of  the 
laser  wavelength  and  on  the  ability  to  measure  it  accurately.  The 
uncertainty,  dx,  of  a displacement  measurement,  x,  due  to  wave- 
length uncertainty,  dX,  is  given  by 


Wavelength  stabilities  on  the  order  of  5 x lO"?  jum  over  an  8-hr  period 
using  an  intracavity  air  spaced  etalon  are  achievable  on  commercially 
available  lasers  (Ref.  14).  However,  changes  in  the  dimensions  of  the 
laser  cavity  and  in  the  index  of  refraction  of  materials  in  the  cavity  due 
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to  temperature  changes  and  mechanical  and  acoustical  vibrations  will 
increase  this  value.  The  temperature  stability  of  an  etalon  is  4 x 10 "8 
nm/°C.  Even  with  these  variations  the  measurement  accuracy  of  the 
interferometer  is  well  within  that  required  to  measure  model  displace- 
ment in  a wind  tunnel. 

A more  important  source  of  error  in  the  displacement  measurement 
is  vibration  of  the  optical  components.  A frequency  shift  induced  in  one 
of  the  heterodyned  light  beams  caused  by  a relative  displacement  in  one 
of  the  optical  elements  will  introduce  extraneous  information  into  the 
displacement  signal.  Therefore,  motion  between  optical  elements  must 
be  kept  to  a minimum.  This  can  be  achieved  with  a compact,  rigid, 
optical  system  mounted  on  a massive  support  isolated  from  acoustical 
and  mechanical  vibrations. 


4.2  DENSITY  EFFECTS 


Since  the  interferometric  technique  measures  the  difference  in 
length  between  two  optical  paths,  the  signal  arising  from  the  mixing  of 
the  two  light  beams  is  a function  of  not  only  the  model  position  but  the 
index  of  refraction  along  each  of  the  two  optical  paths,  L1  and  L2. 

Applying  Fermat's  principle  to  Eq.  (2)  the  relative  variation  of  phase, 
A <t>,  between  optical  paths  is 


A q - - b2 


4 £ 
\ 


(i,t)di 


axt)  di 


(38) 


where  Pj(4,  t)  and  p2(i,t)  are  the  time  and  spatially  dependent  indices  of 
refraction  along  the  optical  paths  Ll  and  L2,  respectively.  Depending 
upon  the  complexity  of  p,  and  the  geometry  of  the  optical  path,  one  can 
either  solve  this  integral  in  closed  form  for  simple  cases  or  a solution 
may  be  obtained  using  numerical  iterative  techniques. 

In  an  operating  wind  tunnel,  density  variations  are  quite  evident. 
The  greatest  variance  is  the  region  between  the  tunnel  wall  and  the 
free-stream  region  where  the  density  is  essentially  constant.  The 
transition  region  .s  called  the  boundary  layer  and  varies  as  the  tunnel 
operating  conditions  vary.  The  density  effect  or  index  of  refraction 
variation  upon  the  interferometric  measurements  must  be  considered. 

In  air,  the  Lorentz-Lorentz  relationship  (p  ^ 1.0)  reduces  to  the 
Gladstone -Dale  relationships  given  by 

b - 1 = Kp  (39) 
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where  p is  the  air  density  and  K is  a proportionality  constant.  Assum- 
ing X = 0.  6328  p,  then  K = 117  x 10-*  ft*  slug-!. 

By  far  the  largest  effects  will  be  due  to  variations  in  the  tunnel  from 
the  time  calibration  measurements  are  made,  through  the  transition  to 
the  particular  tunnel  operating  conditions  where  the  dynamic  data  are  to 
be  taken.  Consider  the  operating  envelope  for  the  AEDC  Propulsion 
Wind  Tunnel  Facility  (PWT)  16  f and  16S  tunnels  as  shown  in  Fig.  18. 


Figure  18.  Total -pressure  operating  envelopes  ior  PWT  Tunnels  16T  and  16S. 

At  the  selected  tunnel  operating  points  the  corresponding  pressures  and 
temperatures  were  used  to  calculate  the  free-stream  density.  Also, 
assuming  the  optical  paths  are  in  the  direction  of  density  gradients,  such 
as  boundary  layers,  Eq.  (38)  reduces  to 

A(t>  ^ ^ At  (40) 

where  At.  is  the  change  in  optical  pathlength  due  to  either  changes  in  ri 
or  in  changes  in  L or  both.  The  apparent  phase  change  cannot  be  sepa- 
rated into  index  of  refraction  changes  and  changes  in  L without  some 
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knowledge  of  the  flow -field  conditions.  In  general,  this  is  easily  obtain- 
ed and  even  a simple  first-order  approximation  will  account  for  a large 
fraction  of  the  variation  in  the  optical  path.  As  an  example,  the  effects 
of  the  optical  path  in  traversing  to  the  center  of  a 16 -ft  test  section 
under  various  tunnel  conditions  are  illustrated.  Given  the  typical  data 
points  in  Table  2,  the  change  in  pathlength  due  to  change  in  index  of 
refraction  from  vacuum  conditions  can  be  calculated  based  on  the  total 
pressure,  typical  operating  total  temperatures,  and  Mach  number.  The 
air  density  in  the  tunnels  can  be  calculated  and  consequently  the  index  of 
refraction  is  obtained.  Using  the  Gladstone -Dale  equation,  and  Eq.  (38), 
the  maximum  variation  in  optical  path  from  vacuum  conditions  was  cal- 
culated. 

Kp  = 77-1  = A Hi  (41) 

Over  an  8-ft  distance,  this  can  be  as  low  as  49  |iim  and  as  high  as  687  pm 
for  extremes  in  the  tunnel  conditions  as  shown  in  Column  1 of  Table  2. 


Corrections  for  variations  can  be  made  in  one  of  two  ways:  (1)  by  an 
on-the-body  retroreflector  to  be  used  as  a reference  distance  or  (2) 
through  calculations  based  on  the  knowledge  of  the  tunnel  conditions.  The 
effects  of  boundary  layers  within  the  tunnel  can  also  be  included  in  the 
analysis.  Semi-empirical  relationships  for  the  density  through  the  bound- 
ary layer  can  be  used  in  conjunction  with  the  Gladstone -Dale  equation  to 
determine  an  effective  index  of  refraction.  An  empirical  relationship  for 
the  density  variation  through  a boundary  layer  is  given  by 


1 . 0.  2 M„2  (1  - r^)2/9 

6o 


(42) 


where  p = the  density  at  position  y from  the  chamber  wall  to  the  free- 
stream  position  60  where  the  density  is  pw.  Assuming  a boundary-layer 
thickness  60  of  1 ft  and  integrating  through  the  entire  boundary  layer,  an 
average  p/p „ can  be  obtained  for  a given  free-stream  Mach  number,  M,,,. 
As  in  the  case  of  free-stream  corrections,  the  variation  over  an  8-ft 
path  was  calculated.  The  maximum  variation  from  free  stream  occurs 
for  M,,,  = 3.0.  This  amounts  to  a 7.4-percent  change  in  previously  calcu- 
lated optical  path  difference  over  8 ft.  By  neglecting  this  correction,  a 
maximum  error  of  13.  7 pm  (87  counts  in  a 1/4  X system)  would  be 
incurred. 
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Density  variations  are  evidenced  in  forms  other  than  in  the  boundary 
layers.  Under  certain  flow  conditions  in  the  transonic  and  supersonic 
operating  region,  shock  waves  are  evidenced.  Associated  with  the  shock 
wave  structure  within  a wind  tunnel  are  strong  density  gradients.  The 
effect  of  these  gradients  upon  light  beams  passing  through  such  a medium 
also  changes  the  optical  pathlength  and  introduces  angular  variations  in 
the  light  path.  In  general,  whenever  the  ray  path  is  not  in  the  same 
direction  as  the  density  gradient,  then  a refractive  bending  of  the  ray  will 
occur.  In  the  case  of  planar  shock  waves,  the  net  result  is  a deviation  in 
the  original  path  for  all  plane  parallel  waves  incident  upon  the  shock  front. 
With  only  a rudimentary  knowledge  of  the  flow  conditions,  the  net  change 
in  optical  path  to  a retroreflec^or  located  upon  a model  can  be  calculated 
with  considerable  accuracy. 

The  effects  of  conical  shocks  upon  incident  radiation  is  similar  to 
wavefront  distortion  by  a lens  with  aberrations.  This  effect  becomes  in- 
creasingly important  with  increasing  cone  angle  and  Mach  number  as  the 
density  variation  across  the  shock  becomes  greater.  Both  the  oblique 
traversing  of  the  optical  path  in  the  shock  system  and  the  focusing  effect 
of  a conical  shock  can  affect  the  integrated  pathlength  to  a model  surface. 

A boundary  layer,  as  a shock  front,  can  also  have  refractive  effects 
on  light  rays  passing  through  it.  The  boundary  layer,  furthermore,  can 
cause  phase  fluctuations  which  are  of  importance  to  interferometric  tech- 
nique. Fluctuating  boundary  layers  can  further  cause  image  dancing, 
scintillation  or  loss  in  image  resolution  which  are  of  particular  importance 
in  imaging  systems  design. 

While  it  is  beyond  the  scope  of  this  report  to  comprehensively  deal 
with  each  of  the  density  related  optical  path  deviations,  it  was  shown  that 
a major  contributor,  i.  e.  , tunnel  condition  changes,  while  significant, 
can  be  calculated  and  used  to  reduce  e.  rors.  Furthermore,  the  density 
effects  can  be  reduced  and  in  some  cases,  eliminated  through  the  use  of  a 
reference  beam  taken  from  the  model  itself.  This  is  due  to  the  matching 
of  optical  path  conditions  accurately  over  the  small  spatial  differences 
along  each  optical  path.  Under  these  conditions,  density  effects  are  mini- 
mized. In  general,  even  when  on-the-body  reference  beams  are  used, 
some  knowledge  of  the  flow  field  will  be  required  to  ensure  accurate  data 
analysis. 
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5.0  SUMMARY  AND  CONCLUSIONS 

An  interferometric  technique  has  been  proposed  for  measuring  model 
displacement  and  q 'formation  in  a wind  tunnel.  A two-beam  noncontact, 
point -measuring  device  with  a resolution  to  one-quarter  of  the  laser  wave- 
length has  been  described.  Data  from  this  device  may  be  recorded  and 
reduced  online  and  are  computer  compatible. 

Relative  measurements  can  be  made  of  the  model  with  the  retro- 
reflectors  mounted  on  the  wing  and  fuselage.  An  absolute  measurement, 
for  model  position  relative  to  a stationary  reference,  can  be  made  by  plac- 
ing one  retroreflector  in  a known  atmosphere.  Retroreflectors  are  used 
to  maintain  optical  alignment  throughout  the  span  of  model  motion. 

A polarization-sensitive  optical  system  reduces  interference  effects 
of  stray  light  from  the  model  and  optical  surfaces  and  separates  the  path 
of  the  illumination  beam  from  the  return  beam.  The  illumination  is  colli- 
mated to  spotlight  areas  on  the  model  containing  a retroreflector  to  effici- 
ently utilize  the  laser  light.  Diameter,  position,  and  angle  of  the  illumina- 
tion beam  are  dictated  by  the  amount  and  direction  of  reflector  motion 
during  the  test.  Optical -to -electrical  signal  conversion  is  performed  by 
an  avalanche  photodiode  which  has  wide  bandwidth,  high  sensitivity,  and 
low  noise.  Because  of  the  efficiency  of  the  photodiode,  the  laser  source 
power  can  be  relatively  low,  thereby  minimizing,  if  not  obviating,  safety 
problems.  The  severe  environment  of  a test  cell  requires  that  the  optical 
system  be  compact  and  rugged  and  isolated  from  both  acoustical  and 
mechanical  noise  and  vibration. 

Displacement  data  can  be  recovered  from  the  15-MHz  frequency- 
modulated  carrier  with  either  a phase  quadrature  digital  pulse  counting 
technique  or  a frequency-to-voltage  conversion  analog  system.  These 
data  can  be  made  available  digitally  for  compute”  analysis  and  as  an 
analog  signal  for  immediate  evaluation  during  and  following  a test. 

The  displacement  interferometer  is  capable  of  resolving  motion  to 
A/4  which  for  a helium-neon  laser  is  158.  2 nanometers.  Maximum  sur- 
face velocity  is  limited  by  the  bandwidth  of  the  electronic  signal  proces- 
^ sor.  A maximum  frequency  deviation  of  5 MHz  from  15  MHz  (10 -MHz 

bandwidth)  limits  the  velocity  to  1.6  meters/sec.  For  oscillatory  motion, 
the  maximum  peak  amplitude -frequency  product  is  25.2  cm-Hz.  Maxi- 
mum displaceme»if’,  on  the  other  hand,  is  limited  electronically  by  the 
number  of  up/down  counters  in  the  digital  data  processing  circuit  and  by 
the  maximum  voltage  and  sensitivity  of  the  analog  system.  With  six  and 
one -half  digits  in  the  readout  of  the  phase  quadrature  digital  system  and 
A/4  resolution,  displacements  of  ±31.64  cm  can  be  measured. 
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ment.  Techniques  are  available  to  compensate  for  these  variations. 
Nevertheless,  future  work  will  include  more  detailed  theoretical  and 
experimental  investigations  into  the  refractive  effects  of  shocks  and 
boundary  layers  on  the  uniformity  of  the  optical  wavefronts  at  the  photo  - 
detector.  Vibration  as  well  as  acoustical  affects  on  the  opiical  package 
will  be  analyzed  to  determine  the  amount  of  isolation  necessary  to  pre- 
vent the  optical  components  from  vibrating  and  introducing  extraneous 
signals  on  the  model  displacement  signal. 

The  two-beam  system  will  be  used  in  actual  tunnel  tests  to  analyze 
the  effects  of  vibration  and  laser  beam  perturbations  due  to  the  tunnel 
flow  conditions.  A multiple-beam  optical  system  will  then  be  designed 
and  built  to  increase  the  capability  for  simultaneous  measurements  of 
model  attitude  and  deformation. 

Electronic  signal  processors  for  the  multiheam  system  will  be 
fabricated  with  improvements  incorporated  in  bandwidth  and  response 
as  required  for  specific  applications.  The  use  of  a microcomputer  for 
the  online  acquisition  and  display  of  data  from  a set  of  multibeam  read- 
outs will  be  studied. 

The  possible  replacement  of  the  retroreflectors  with  highly  reflect- 
ing tape  or  a diffuse  surface  will  be  analyzed.  This  would  make  model 
preparation  much  simpler.  Me?r°  of  reducing  the  number  of  optical 
components  will  also  be  studied.  This  includes  scanning  single  or 
multiple  laser  beams  across  the  surface  and  collecting  the  backscatter- 
ed  light  by  a single  lens.  Frequency  multiplexing  of  the  signals  will  be 
studied  as  a possible  means  of  reducing  the  number  of  photodetectors  in 
the  multiple -beam  system. 
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V = 1,2. 
m 

A(0.) 
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di»  ds 


Eo’  ER 
Af 

fl*  f2 

fc 

f(t) 

Gn  = 1,  2.  . . 
g(t) 


NOMENCLATURE 

Zero-crossing  detected  signal  of  Pj 
Electrical  signal  amplitude 
Effective  aperture  of  retroreflector 
Delayed  A waveform 
Zero- crossing  detected  signal  of  P2 
Delayed  C waveform 
Down  count  pulse 

Beam  diameter  leaving  optics  package 

Diameter  of  retroreflector 

Beam  diameter  prior  to  beam  expanding 
optics 

Diameter  of  returning  beam  inside  optics 
package 

Signal  from  34-MHz  and  49-MHz  oscillator, 
respectively 

Amplitude  of  object  and  reference  beam 
Bandwidth  of  electronics 
Focal  lengths 

Carrier  frequency  (15  MHz) 

Signal  frequency 

Gain  of  electronic  components 

Sum  and  difference  electronic  signal 
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Distance  between  output  optics  and  model 
surface 


*0'  *R 


Ll*  L2 


M^t),  M2(t) 


NEP 


P^t),  P2(t) 


r(t),  r^t),  r2(t) 


Intensity  incident  on  retroreflector 

Intensity  of  object  and  reference  beam 

Photodetector  output 

Proportionality  constant 

Lateral  motion  of  model 

Optical  path  length  of  object  beam  and 
reference  beam 

Lateral  motion 

Signal  after  mixing  with  reference  oscillator 
signal 

Free- stream  Mach  number 

Number  of  light  beams  generated 

Noise  term  due  to  detector  and  optical 
noise 

Noise  equivalent  power 

Number  of  frequency  division  stages  in 
analog  processor 

Lowpass  filtering  of  M^(t)  and  M2(t) 

Distance  from  point  of  rotation  to  retro- 
reflector 

Reference  signals  (see  Fig.  7) 
Retroreflector  spacing 
Frequency -modulated  signal 
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U 

Av(t) 

x(t) 

y 

a 

3 


7 
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Beam  spacing  prior  to  beam -expanding 
optics 

Up  count  pulse  train 

Velocity  difference  between  two  retro - 
reflectors 

Reflector  motion 

Distance  from  chamber  wall  (inside  of 
boundary  layer) 

Model  rotation  angle 

Air-glass  interface  loss 

Retroreflector  efficiency 

Polarization  rotation  angle  upon  reflection 

Boundary -layer  thickness 

Bragg  cell  light  deflection  efficiency 

Index  of  refraction 

Angle  of  incidence  of  light  beam 

Maximum  model  illumination  angle 

Optical  wavelength 

Density 

Phase  term  of  object  and  reference  beam 

Optical  carrier  frequency  or  (l>0  - u>R) 

Angular  frequency,  object  beam 

Angular  frequency,  reference  beam 
Frequency  deviation 
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